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Abstract 

Recent experimental studies on the instability appearance of capillary jets have revealed the ca- 
pabilities of linear spatiotemporal instability analysis to predict the parametrical map where steady 
^H ! jetting or dripping takes place. In this work, we present an extensive analytical, numerical and ex- 

O ' perimental analysis of confined capillary jets extending previous studies. We propose an extended, 

^.^' accurate analytic model in the limit of low Reynolds flows, and introduce a numerical scheme to 

^. 
CL|, predict the system response when the liquid inertia is not negligible. Theoretical predictions show 
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a remarkable accuracy with results from the extensive experimental exploration provided. 



I. INTRODUCTION 

An important number of advances in the understanding of the physics of the controlled 
formation of biphasic flows {i.e. emulsions, aerosols, microbubble dispersions, wet foams, 
etc.) have come after the incorporation of spatiotemporal analysis in the dynamics of cap- 
illary jets. A general successful paradigm in such controlled formation mentioned is the 
co-flow of the two phases involved, where a first fluid stream from a given source, usually a 
capillary tube, is surrounded by a second immiscible fluid stream. Basically, the first fluid 
stream may exhibit two possible behaviors: (i) it gives up blobs (droplets or bubbles) of 
sizes comparable to the typical source dimensions (usually the tube outlet diameter), or (ii) 
it yields a continuous stream in the form of a capillary jet whose diameter becomes inde- 
pendent of the typical source dimensions. These system modes are usually referred to as 
"dripping" or "jetting", respectively [l|, |2|,|3|, y,|5|. Id, LZf. In this later case, after the capillary 
jet breaks up into droplets a distance of several diameters downstream of the source, the 
eventual output is a stream of droplets or bubbles whose size is commensurate with the jet 



diameter. When subject to flow under certain geometrica 



co nfig urations (like flow focusing, 



P, y, [lO| ) or electrohydrodynamic actions (electrospray, ll|, [l2]), the first fluid stream can 
be reduced to the form of a steady jet of very small diameter compared to the typical size 
of the source. In general, producing a biphasic flow based on jetting offers very significant 
advantages over dripping because (i) a fixed geometry can be tuned to yield a wide range of 
selectable droplet sizes based on relative flow rates only, (ii) the resulting droplet sizes can be 
sometimes extremely small compared to the given geometrical dimensions of the processing 
device, and (iii) under certain parametrical circumstances, the droplet size distribution can 
exhibit a low size dispersion. Thus, for each device geometry, mapping the different behav- 
iors mentioned in the parametrical space of operation becomes a fundamental engineering 
endeavor. 

An extremely simple yet successful device to produce such controlled biphasic flows is a 
concentric capillary tube arrangement [13l. Il4l |. Using such a geometry, a smaller diameter 
tube coaxially positioned inside of another wider confining tube releases a stream of a first 
fluid, usually a liquid, while a second fluid flows concentrically in the same axial direction 
(for example, see Fig. [T]). This axisymmetric arrangement is rather common to many mi- 



crofluidic systems. Very recently, Guillot et al. lj| have contributed very significantly to the 



physical knowledge of such a simple yet fundamental system, presenting a combined exper- 
imental and theoretical study with a remarkable degree of agreement at many parametrical 



regions, in spite that their theoretical model proposed assumed important simplifications: 
both the basic flow and its perturbations were averaged in the radial direction for both 
streams, and inertia was neglected through the whole study (justiflable in microfluidics) . 
These simplifications naturally become shortcomings in some cases where predictions devi- 
ated from experiments. In this work, we aim to investigate whether an axisymmetric model 
which does not make use of these simplifications improves its predictive power over the sim- 
plified one. To do so, we split our study into two subsequent degrees of approximation to the 
real physics: firstly, we assume the perturbed fiow three dimensional axisymmetric under 
negligible inertia, and secondly, we add inertia. The first study still admits analytical treat- 
ment in spite of the extremely complexity of the resulting mathematical expressions, which 
should be handled by a computer assisted algebraic manipulation {e.g. Mathematica®) , 
and the second study is numerically tackled. Our theoretical results are compared to an 
extensive collection of experiments, an important part of which is here presented for the 
first time. As expected, the theoretical agreement with experiments improves significantly. 
Still, some remaining disagreements are discussed in detail, possibly shedding some light on 
the dynamics of these fiows for future studies. 

II. MATHEMATICAL FORMULATION 

We studied the spatial response to small perturbations of an infinite cylindrical capil- 
lary liquid jet of radius Ri confined in a cylindrical coaxial pipe of radius R2 > Ri- The 
liquid jet of density pi and viscosity /ii is surrounded by another liquid of density p2 and 
/i2- The governing equations are the incompressible Navier-Stokes equations in cylindrical 
coordinates, {z,r,9), which are made dimensionless by using the radius of the jet, Ri, the 
velocity at the liquid-liquid interface, U, and inertia piU^, as characteristic quantities. This 
problem admit an exact steady unidirectional solution which only depends on two dimen- 
sionless parameters; the viscosity ratio fi = fi2/pi and the quotient between the pipe wall 
radius an the radius of the liquid jet, R = R2/R1. The velocity field {U,V) of this basic 
solution is given by: 

f/i(r) = 1+^1^, Vi{r) = (0<r<l), 

(1) 
U2{r) = 1 + ;^^, V2{r) = {l<r<R), 

where subscript 1 (2) denote the inner (outer) liquid velocity fields. To study the stability 
of the fiow, we will introduce a small unsteady axisymmetric perturbation in the velocity 



and pressure fields of tlie form, 



{u'n,Vn,p'n}{r,z) = {un,Vn,Pn}{r)exp[i{kz -ujt)], (n=l,2), (2) 

where u is the dimensionless wave frequency, k is the dimensionless wave number, u and v 
are the axial and radial velocities respectively, and p is the pressure field. In addition to 
this, the radial position of the liquid-liquid interface is also perturbed in the form 



ri{z,t) = 1 + CoCxp [i{kz — ut)] 



(3) 



where Co is the small amplitude of the perturbated interface ( Co « 1). 

Substituting the basic flow given by ([T]) and the perturbed one given by ([2]) and Q into 
the Navier-Stokes equations and after linearizing we get the following set of equations: 
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where 5„2 is the Kronecker delta. In the above equations. Re represents the Reynolds number 
of the inner liquid, defined as Re = piURi/pi, and p = P2/P1 is the density ratio. Equations 
(jl])-([7j) must be solved subject to the following boundary conditions: 



At the axis, r=0, regularity conditions: 



t^i = 0. 



dui 
dr 



0. 



At the liquid-liquid interface, r = 1, continuity of both the velocity filed and tangential 
stresses, and normal stresses balance with capillary forces which yield: 



Vi = V2, Ui 



dui 



duo 



f^- 
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(9) 



Re dr Re dr We 
being. We = piU'^Ri/a, the Weber number based on the inner liquid density, the 
radius and the velocity of the jet at the interface and the liquid-liquid surface tension 



At pipe wall, r = R, no slip conditions are applied: 



V2 = U2= 0. 



(10) 



Finally, by making use of the linearized kinematic condition for the interface we get: 



iuCo — ik + vAr = \] 



(11) 



Here, we will carry out a spatiotemporal stability analysis, which entails assuming both 
the frequency u and the wave number k = kr + iki as complex. Using a well established 
criterion |14j. Iisl . Ilq . we say that a mode is convectively unstable if its spatial growth rate, 
—ki, and its group velocity, Cg = du/dkr, are both positive for some frequency range. 



A. Creeping flow limit, Re — > 0. 



System of equations (H1)-( ITT1) can be substantially simplified when viscous effect are dom- 
inant {Re — > 0). In this limit the stability problem is governed by the following equations: 

(12) 
(13) 
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Here, p = Re ■ p is the re-scaled pressure amplitude. 



At the axis, r=0: 



^1 = 0, 



dui 
dr 



(16) 



At the liquid-liquid interface, r = 1: 
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dr 
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(17) 



where now, Ca = We/ Re, is a capillary number. 



• At pipe wall, r = R: 

V2 = U2= 0. (18) 

The linearized kinematic equation for the interface reads: 



iujCo - ik + vi{r = 1) = 0. (19) 

Observe that while the complete problem is governed by five dimensionless parameters 
p, /i, Re, We and R, the creeping flow limit is governed by just three, //, Ca and R. 

B. Numerical scheme for the complete system 

It proves convenient to rewrite Equations fHl)- flTTl) in the form 

[Li + A;L2 + — L3]S = 0, (20) 

JrCe 

where Li, L2, L3 and L3 are complex matrices which depend on r, p, /i and u and where 

S = [ui{r),Vi{r),pi{r),U2{r),V2{r),p2{r),Co]. To solve numerically this equation we have 

used a Chebyshev spectral collocation technique based on that developed by Khorramin 



17|, for the stability analysis of swirling flows in pipes. This code have been successfully 



used in the past to analysis the stability of low density and viscosity fluid jets and spouts 
in unbounded coflowing liquids [l8|. To implement the spectral numerical method, (I2U]1 
is discretized by expanding S in terms of truncated Chebyshev series in each sub-domain. 
The interval < r < 1 is discretized and mapped into the Chebyshev polynomial domain 
— 1 < ^ < 1 using the algebraic transformation 



r 



_ (1 - 0) 
' 2 



(j = l,-,iVi) (21) 



r, = l + (i?-l)^^^, (j = l,..,iV2) (22) 



while the interval 1 < r < i? is mapped into — 1 < ■C ^ 1 using 

2 

where A^i and N2 denote the numbers of collocations points used in the radial discretization of 
;he two domains. The discretized system is solved using the linear companion matrix method 
19|. The resulting linear eigenvalue problem is solved using a eigenvalue solver subroutine 
DGVCCG from the IMSL library, which provides the entire spectrum of eigenvalues (k) 
and eigenf unctions (S). Spurious eigenvalues were ruled out by comparing the computed 
spectrums obtained for different values of the numbers A^^i and N2 of collocation points. 
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C. Analytical dispersion relation for the creeping flow limit 

Equations (I12til4l) can be solved by the use of the Stokes stream function \l/n such that 

Un = --r ir'^'^n) , Vn = -tkr^i/r, (23) 

r ar ^ ' 

The solutions of p„ and ^„ verifying fll2m4p are: 

pi(r) = -2iA;ai/o(A:r), (24) 

P2(r) = -2zA;/i[6i/o(A;r) + 62A'o(A;r)], (25) 

*i(r) = AiJi(A;r)/r + a^h{kr), (26) 

*2(r) = [5iJi(A;r) + B^K^{kr)\ /r + hih{kr) + b2Ko{kr). (27) 

The seven conditions (11611191) resolve the six factors {ai,bi,b2, Ai, Bi, B2} and provide 
the dispersion relation: 

where Di = Di{k, /i, R) and A'^i = Ni{k, fi, R). The lengthy resulting expressions, not given 
here for simplicity, have been calculated and the problem handled using Mathematica®. 

III. EXPERIMENTAL 

To obtain a jet in a cylindrical glass capillary of inner radius R2 we use as a nozzle a 
glass capillary of square cross-section with a tapered end (see Fig. [1]). The outer diagonal of 
this square capillary is very close to the inner diameter of the cylindrical tube which ensures 



good alignment and centering of the nozzle in the cylindrical capillary [13|]. In this study 
R2 is equal to 275 or 430 /im, whereas the radius of the tapered orifice of the square tube 
is set between 20 and 50 /xm using a pipette-puller and microforge set up. Syringe pumps 
are used to inject the fluids. Inner fluid of viscosity /ii is injected at rate Qi in the square 
capillary and the outer fluid of viscosity /i2 is injected at a rate Q2 through the cylindrical 
capillary. This leads to a coaxial injection of both fluids at the tapered orifice. 

Flow patterns are observed through an optical light microscope and pictures are recorded 
with a fast camera. They vary significantly with operational (Qi, Q2), geometrical (-R2), 
and system parameters (/ii, yU.2, a). Fluids used in this study are water at 1 mPa-s, mixtures 
of water and glycerine at 55 and 650 mPa-s, silicone oil at 235 niPa-s and hexadecane at 



FIG. 1: Microfluidic device used to perform the experiments. In inner fluid is injected at rate Qi 
through the tip of a tapered square capihary inside a cyhndrical capihary in which an outer fluid 
is injected at rate Q2- 
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FIG. 2: Map of the flow behavior in the {Qi, Q2) plane. The droplet regime comprises droplets 
smaller than the capillary (o) and plug- like droplets confined by this capillary (•). Jets are observed 
in various forms: jets with visible peristaltic modulations convected downstream (D), wide straight 
jets (■) that are stable throughout the ~ 5 cm long channel, and thin jets breaking into droplets 
at a well defined location (0). Parameters are: R2 = 275 /um, inner viscosity Hi = 55 cP, outer 
viscosity fi^ = 235 cP, surface tension a = 16 mN/m. 

3 mPa-s. The surface tension between the fluids are changed by adding sodium dodecyl 
sulfate (SDS) in aqueous solutions and Span 80 in hexadecane. 

Figure [2] displays the typical outcome of an experiment where the flow rates are varied for 
a given system (here the inner solution is a 50 % in weight glycerine in water solution with 
SDS for which /ii = 55 mPa-s and the outer one is a silicone oil for which /i2 = 235 mPa-s). 
A droplet regime is found for low Qi, with either droplets emitted periodically right at the 
nozzle - symbol (o) - or non spherical plug-like droplets resulting from the instability of 
an emerging oscillating jet (•). Jets are found in the bottom right corner of Fig. [2] with 
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different visual aspects: wavy jets with features that are convected downstream (D) and, for 
larger values of Qi, straight jets (■) that persist throughout the cylindrical capillary. For 
large values of the external flow rate Q21 we observe what we call jetting: thin and rather 
straight jets (0) that extend over some distance in the capillary tube before breaking into 
droplets at a well-defined and reproducible location. This "jet length" increases with Qi for 
a fixed value of Q2- Similar "dynamic phase diagrams" have been reported for similar and 
more complex microfiuidic geometries 9|, |20|, l2l|, |22|, |23|, |2J, |25 1 . 



IV. RESULTS 

In order to compare the numerical results with the experimental ones we have first to 
relate the dimensional parameters used in the experimental setup [i.e. the imposed liquid 
flow rates Qi and Q2 of the inner and outer liquids, respectively, i?2, cr and the densities 
and viscosities of both liquids) with our dimensionless parameters i?e. We, R, /i and p. To 
do that we need to obtain the radius of the jet, -Ri, and the velocity at the interface, U, 
from the imposed flow rates and geometrical constrain. By assuming that the flow is fully 
developed with a velocity profile given by ([T]) we get: 

R, = R2{l + Q[l + il + p/QY^']}-'^\ U=^^^^y (29) 

where Q = Q2/Qi- From these expressions, for a given pair of liquids and a certain capillary, 
one can obtain the theoretical mapping in a {Qi, Q2} plane of the absolutely and convectively 
unstable regimes, corresponding to dripping and jetting regimes, respectively. 

Figures ([3|)_provide a comparison of (i) the approximate analytical model proposed by 
Guillot et al.\lj^, (ii) the exact analytical model here proposed for Re = 0, and (iii) the com- 
plete model including inertia, here numerically solved using Chebyshev collocation methods. 

V. DISCUSSION 

1. Low outer flow rate values 

This region deals with the transition from the plug regime to the wavy/straight jet regime. 
Jets obtained in this region are highly confined by the walls. In our experimental conditions, 
Reynolds number is always lower than one in this region. 
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FIG. 3: Map of the flow behavior in the {Qi, Q2) plane for various geometries and fluids. Gray 
symbols correspond to droplets regimes and black symbols to jets regimes. Lines are the transition 
predict by the model without adjustable parameters. The red line corresponds to the linear analysis 
perform using the averaged (lubrication) approximation (Guillot et al. 2007), the blue line to the 
exact analytical model without inertia and the blue line to the numerical one taking into account the 
inertia. Comparisons between experimental and theoretical transitions are in good agreement for 
all the experimental conditions presented here, (al) is obtained for R2 = 275 ;um, ^li = 55 mPa-s, 
/i2 = 235 mPa-s and cr = 24 mN/m\ (a2) for R2 = 430 /iin, /ii = 55 mPa-s, ^2 = 235 mPa-s 
and (T = 24 mN/m; (a3) is o btained with R2 = 275 /im, /ii = 55 mPa-s, 112 = 235 mPa-s and 
o" = 16 mN/m] (bl) is obtained with R2 = 275 /.im, ^i = \ mPa-s, /i2 = 3 mPa-s and cr = 12 mN/m\ 
(b2) is obtained with R2 = 275 /xm, ^ui = 1 mPa-s, /i2 = 3 mPa-s and a = 0.12 mN/m] (cl) is 
obtained with R2 = 275 ;um, /.ii = 235 mPa-s, ^2 = 55 mPa-s and cr = 24 mN/iri] (dl) is obtained 
with R2 = 275 ;Um, /ii = 650 mPa-s, ^2 = 235 mPa-s and cr = 24 mN/m. 

In all the experiments, the three models predict almost the same transition at low values 
of Q2 with a very good degree of agreement. As expected, inertia does not play a major 
role for this range of flow rate, and the differences between the complete model with or 
without inertia cannot be noticed. However, the analysis obtained in the framework of the 
lubrication approximation gives a transition at slightly lower (higher) values of Qi when 
/i >1 (yU <1). Results obtained through this approximate analysis seems to be a little bit 
more accurate that the ones obtained with the whole model. The main difference between 
the lubrication approximation used to perform the linear analysis [M] , which holds for high 
conflnement at any viscosity ratio, and the "exact" models is the former does not consider the 
normal viscous stresses balance at the interface, while the later ones do so. Interestingly, at 
the time of the experimental comparison, perturbation amplitudes may become comparable 
to the thickness of the outer layer at high conflnement, even though they remained small 
compared to the jet radius; in these cases, the accuracy of exact linear models is overcome 
by the approximate averaged model. Informally speaking, the averaged -more "forgiving" - 
model seems to yield a better approximation for visible perturbation amplitudes, while the 
"exact" ones -more sensitive- should depart more signiflcantly from experiments owing to 
non- linear effects for high conflnement (flnite perturbation amplitudes). 

Thus, at low Reynolds number and for high value of conflnement, the various analyses 
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presented in this paper allow to predict quantitatively the transition between a droplet and 
jet regime. In spite of the number of approximations, the averaged model appears to be a 
simple and powerful tool to predict the transition in this range of low outer flow rate. 

2. High outer flow rate values 

This range of flow rate concerns the region where, by changing the flow rate, we can 
obtain a transition between a droplet and a jetting regime. Jets obtained in this region are 
thin and hardly affected by the walls. Note that the Reynolds number can be higher than 
the unity in some of our experimental conditions. 

When we compare the experimental observations with the results obtained by the different 
models it clearly appears that the complete "exact" models provide a better prediction of 
the transition between the droplet and the jetting regime than the averaged model. We can 
also notice the effect of inertia in Fig. [3](b2) where the two lines clearly split at high values 
of Q2- Unfortunately, we do not have enough experimental data to validate the fact that 
inertia promotes transition to dripping for smaller outer flow rates. 

The averaged model fails for low degree of confinement. This is mainly due to the fact 
that the radial component of the velocity can not be neglected over the axial velocity and 
that the recirculations inside the perturbed jet have to be taken into account, as well as 
the normal stresses at the interface. In other words, the averaged lubrication approximation 
used is not valid in this case. The range of R in which the approximation still holds depends 
on the viscosity ratio between both fluids. 

When the transition towards jetting is studied, the complete model clearly gives better 
results and allows to predict quantitatively the transition between dripping and jetting. 
The effect of inertia will have to be checked experimentally to confirm the prediction of Fig. 
Hbl). 

VI. CONCLUSIONS 

In this work, we have studied the stability of confined capillary jet at low Reynolds 
numbers by considering the transition from a steady jet to a dripping system in terms of 
a convective to absolute instability transition. We propose both analytical and numerical 
analyses to describe extensive experimental results obtained in a cylindrical geometry. Pre- 
dictions obtained in high degree of jet confinement are in very good agreement with all the 

12 



analysis presented in this paper. The hnear analysis proposed in [ij] using the lubrication 
approximation seems to be a powerful and rather simple tool to predict the transition in 
this range of confinement. However this analysis fails for low degree of confinement and, in 
this latter case, it is necessary to go through a complete analysis of the flow to get a good 
agreement between the experimental results and the predictions. The effect of inertia has 
also been introduced in this work. Inertia seems to promote droplet formation at lower outer 
flow rates, but we do not have enough to data to confirm experimentally this prediction. 
Thus, further studies should be performed in this range of flow rate to confirm this effect. 
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